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Unitarization schemes can be reduced to non-linear equations which saturate at small b the evo- 
lution of the elastic amplitude with s, and which mimic parton saturation in the non-pcrturbative 
regime. These equations enable us to study the effect of saturation on total and elastic cross 
sections for various models, and to evaluate the uncertainties on atot and p(s,t) at the LHC in 
the presence of a hard pomeron. 
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1. Introduction 

The most important results on the energy 
dependence of diffractive hadronic scattering 
were obtained from first principles (analytic- 
ity, unitarity and Lorentz invariance), which 
lead to specific analytic forms for the scatter- 
ing amplitude as a function of its kinematical 
parameters - s, t, and u. One of the impor- 
tant theorems is the Froissart-Martin bound 
which states that the high-energy cross sec- 
tion for the scattering of hadrons is bounded 
by otot ~ log 2 (s/so), where sq is a scale fac- 
tor. 

Experimental data reveal that total cross 
sections grow with energy. This means that 
the leading contribution in the high-energy 
limit is given by the rightmost singularity, 
the pomeron, with intercept exceeding unity. 
In the framework of perturbative QCD, the 
intercept is also expected to exceed unity by 
an amount proportional to a s At leading- 
log s, one obtains the leading singularity at 
J — 1 = 12 log2(a s /7r). In this case, the 
Froissart-Martin bound is soon violated. 

In a recent study, 1,2 we have found 
that forward data, including total cross sec- 



tions and the ratios of the real part to the 
imaginary part of the scattering amplitude, 
are well fitted by a combination of a soft 
pomeron (which would be purely non pertur- 
bative) and a hard pomeron. The inclusion 
of these two pomcrons, together with the use 
of the integral dispersion relations, 3 and the 
addition of the sub-leading meson trajecto- 
ries, leads to a successful description of all pp, 
pp, 7r ± p, K ± p, 7p and 77 data for y/s < 100 
GeV. Indeed, the fast growing hard pomeron 
leads to a violation of unitarity for values of 
y/s around 1 TeV. 

Now the saturation processes are ac- 
tively studied in high-energy physics. They 
are connected with the saturation of the 
gluon density at small x, which may be re- 
lated to the Black Disk Limit (BDL). How- 
ever, one should note that the process of 
saturation has mostly been examined in the 
framework of the dipole model. 4 But, as 
noted by the Golec-Biernat and Wusthoff, 5 
saturation in that case suppresses the soft 
contributions and the hard Pomeron plays 
the main role. That is why the dipole model 
describes so successfuly hard processes (such 
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as vector-meson production). Contrarily, 
saturation in the BDL sense tames the hard 
Pomeron and increases the role of the long- 
distance processes. It is the latter that we 
shall study in this contribution. 

The unitarity condition can be imple- 
mented via several different prescriptions. 
Two of them arc based on particular solu- 
tions of the unitarity equation. The scatter- 
ing amplitude in the impact parameter rep- 
resentation is defined as 

/>oo 

T(s,t) = i bdbJ (bA)f(b,s). (1) 
Jo 

with f(b,s) < 1. It satisfies the unitarity 
equation 

Imf(s,b) = [Im(f(s,b)} 2 + [Ref( S ,b)} 2 

+ Vinel(s,b). (2) 

One of the possibilities is obtained in the 
U— matrix approach 6 ' 7 : 



f(s,b) 



U(s,b) 



(3) 



1 - i U{s,b) 

The second possible solution of the uni- 
tarity condition corresponds to the eikonal 
representation 

f(s,b) = [l-exp(-xM)]. (4) 

One often takes the eikonal phase in fac- 
toriscd form 



X (s,b) = h(s) f(b), 



(5) 



and one supposes that, despite the fact that 
the energy dependence of h(s) can be a power 
h(s) <~ s A , the total cross section will satisfy 
the Froissart bound a to t < a log 2 (s) We 
find in fact that the energy dependence of the 
imaginary part of the amplitude and hence 
of the total cross section will depend on the 
form of f(b), i.e. on the t dependence of the 
slope of the elastic scattering amplitude. If 
f(b) decreases as a power of 6, the Froissart- 
Martin bound will always be violated. In the 
case of other forms of the b dependence, a 
special analysis 8 is required. 



2. Non-linear equations 

The problem of the implementation of uni- 
tarity via saturation is that the matching 
procedure seems arbitrary. Hence we consid- 
ered a different approach to saturate the am- 
plitude. It is connected with the non-linear 
saturation processes which have been consid- 
ered in a perturbative QCD context. 9 ' 10 Such 
processes lead to an infinite set of coupled 
evolution equations in energy for the corre- 
lation functions of multiple Wilson lines. 1 
In the approximation where the correlation 
functions for more than two Wilson lines fac- 
torise, the problem reduces to the non-linear 
Balitsky-Kovchegov (BK) equation. 11 ' 12 

It is unclear how to extend these results 
to the non-perturbative region, but one will 
probably obtain a similar equation. In fact 
we found simple differential equations that 
reproduce either the U-matrix or the eikonal 
representation. 

We shall consider saturation equations of 
the general form 13 

dN/d£ = S(N) (6) 

We shall impose the unitarity conditions 
TV — > 1 as s — > oo and dN/dt; — > as 
s — > oo. We shall also assume that S(N) has 
a Taylor expansion in TV, and that the first 
term only the hard pomeron Nb are = f(b)s A , 
where A is the pomeron intercept minus 1. 
Similarly, we fix the integration constant by 
demanding that the first term of the expan- 
sion in s A reduces to Nbare ■ We then need to 
take S(N) = AN + 0(N 2 ). The conditions 
at s -> oo then give S(N) = A(N — N 2 ) as 
the simplest saturating function. The result- 
ing equation 

dN/dy = A(N - N 2 ), (7) 

where y — log(s) ~ log(l/x) This equation 
has the solution 

_ f(b)s A 

N -m^T~i (8) 

One may then wonder whether other uni- 
tarization schemes are possible. We give here 



February 2, 2008 7:30 WSPC/Trim Size: lOin x 7in for Proceedings 



ihcp06d 



a general algorithm to build such schemes: 
dN dN bari 



dy 



dy 



-(l-N) 



(9) 



with Nb are the unsaturated amplitude. This 
will trivially obey the conditions at s — > oo, 
and saturate at N = 1. It has as a solution 

N(b,s) = 1 - exp(-AWe(M)) (10) 

This is equivalent to the eikonal representa- 
tion with N bare (b, s) corresponding the Born 
aplitude in the impact representation. The 
corresponding equation can also be written 
in the form 

dN 

— = A N bare exp[-N bare (b,s)]. (11) 
dy 

One can also find new unitarization schemes, 
e.g. from the equation 

dN / dy = (i _ N ^ (12) 



dy 



one gets 



N(b,s)=t&nh[N bare (b,s)}. (13) 

On Fig.l, the result of saturation is shown 
for these three unitarization procedures. It 
is clear that despite the essentially difference 
in the forms of the non-linear equations, the 
process of the saturation is very similar. 




8 10 12 



Fig. 1. The saturation of N(y) for different non- 
linear equations (hard line - eq.(ll, long- dashed line 
- eq. (7), and short-dashed - eq. (12)) 



3. Minimum saturation of the 
Soft plus Hard Pomeron 
model 

Let us take, as an example, the Soft plus 
Hard pomeron model 2 ' which includes two 
simple poles to describe pp and pp scattering 
at high s. In this case, the pp-elastic scatter- 
ing amplitude is proportional to the hadron 
form factors and can be approximated at 
small t by: 

T(s,t) ~ [ Jn (s/s ) ei e a 'i * lo s( s /«o) 
+h 2 (s/s y 2 e a * * lo s( s /«o)] F 2( t y ( 14 ) 

where h\ = 4.7 and h 2 = 0.005 are the cou- 
pling of the soft and hard pomerons, and 
ei = 0.072, a[ = 0.25 GeV~ 2 , and e 2 = 0.45, 
a' 2 = 0.20 GcV~ 2 are the intercepts and the 
slopes of the two pomeron trajectories, and 
where s contains implicitly the phase factor 
exp(— in/2). F 2 (t) is the square of the Dirac 
clastic form factor. It can be approximated 
by sum of three exponentials. 

We then obtain in the impact parameter 
representation a specific form for the profile 
function T(b,s). 8 One finds that for y/s w 
1.5 TeV and at small b, T(b,s) reaches the 
black disk limit. 

Note that one cannot simply cut the 
profile function sharply as this would lead 
to a non-analytic amplitude, and to specific 
diffractive patterns in the total cross section 
and in the slope of the differential cross sec- 
tions. Furthermore, we have to match at 
large impact parameter the behaviour of the 
unsaturated profile function. We have to use 
some specific matching patters which softly 
interpolate between both regimes. The in- 
terpolating function gives unity in the large 
impact parameter region and forces the pro- 
file function to approach the saturation scale 
b s as a Gaussian. 

We show in Fig. 2 the possible be- 
haviours of the total cross section at very 
high energies, depending one = a(0) — 1 and 
on the unitarization scheme. We also show 
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there the result of a simple cikonalisation, 
where we took 

T pp {s, t = 0) ~ J d 2 b[l - e h " kG(s ^} (15) 

and where h ei k — 1.2 was chosen so that 
the values cr tot determined by the overlapping 
function (14), and by cikonalisation (15) are 
equal at y/s = 50 GeV. 
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Fig. 2. The total cross section at ^(s) = 1.8 TeV 
and \/s = 14 TeV (lower two and upper two lines) 
for the saturated amplitude (plain and short dashes) 
and for an eikonalized amplitude (long dashes and 
dash-dots). 



0.3 































X 





0.00 0.05 0.10 0.15 0.20 0.25 



-t (OeVJ) 

Fig. 3. The ratio of the real to the imaginary part 
of the amplitude as a function of t for the saturated 
amplitudes at various energies: 100, 500, GeV, and 
5, 14 TeV (plain, short-dashed, dot-dashed lines cor- 
respondingly 

To demonstrate the presence of satura- 
tion at the LHC, one can also consider the 
slope of the differential elastic cross section, 
which we show in Fig. 3. We see that satura- 
tion increases its value at small t (this is in 
fact unavoidable for any saturation scheme) , 
and predicts a fast drop around \t\ — 0.25 
GeV 2 , when one enters the dip region. 

In conclusion, we have shown that the 
most usual unitarization schemes could be 



cast into differential equations which are 
reminiscent of the saturation equation . n > 12 
Such an approach can be used to build new 
unitarization schemes and may also shed 
some light on the physical processes respon- 
sible for the saturation regime. 
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